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Abstract 

The hadron/pion ratio is calculated in 200 GeV AuAu collisions at midrapid- 
ity, applying pQCD and non-universal transverse-momentum broadening. Ar- 
guments are presented for such non-universality, and the idea is implemented 
in a model, which reproduces the main features of the centrality dependence 
of the hadron/pion ratio in AuAu collisions. The model also reasonably de- 
scribes the qualitative difference between the recently-measured dAu nuclear 
enhancement factors for pions and charged hadrons. 
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Nuclear collisions are studied at unprecedented energies at the Relativistic Heavy Ion 
Collider (RHIC), revealing several new phenomena embedded in a large amount of high- 
quality data, most aspects of which are consistent with general expectations [1]. A lot of 
attention is centered on final-state particles (secondaries) with high transverse momenta 
(pt), where perturbative quantum chromodynamics (pQCD) should have good predictive 
power. In this regime, a suppression of total charged hadron and pion production has 
been observed relative to a nucleon-nucleon (NN) reference in AuAu collisions [2,3]. The 
suppression can be described by pQCD calculations after incorporating final-state partonic 
energy loss in dense matter [4-6] , or using a model of initial-state gluon saturation [7] . The 
most recent dAu data [8] provide important information for understanding RHIC physics. 

Furthermore, the PHENIX collaboration reports an anomalous enhancement of the high- 
Pt proton-to-pion (p/tt) ratio in AuAu collisions at ^/s — 130 GeV [9] and = 200 GeV 
[10]. At 200 GeV, PHENIX also measured the ratio of charged hadrons = {h+ + /i")/2) 
to neutral pions, and found it enhanced in central collisions, presumably due to the enhance- 
ment of proton production. Different versions of coalescence (recombination) of partons from 
the Qurak-Gluon Plasma were proposed to understand the enhanced p/n ratios [11]. Nev- 
ertheless, the measured azimuthal correlations for high-p^ charged particles [12] and the 
binary-collision scaling of proton production gleaned from comparing different centralities 
[10] indicate that the production from hard scattering remains important. Enhanced pro- 
ton/pion and kaon/pion ratios have also been observed in lower-energy pA collisions [13]. 
Therefore, other effects could also contribute to the enhanced particle ratio at RHIC. In this 
paper, we will try to provide an alternative explanation for the particle ratio data at RHIC 
based on pQCD. 

Another experimentally favored quantity, where differences between and tt secon- 
daries manifest themselves, is the "nuclear modification factor", Rab- This ratio is designed 
to display the effects arising in an actual AB nuclear collision, relative to a hypothetical 
collection of independent NN collisions. Earlier experimental information on Rab for AuAu 
has recently been augmented by data obtained in the dAu run to provide a crucial reference 
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[8], and the nuclear modification factor was also addressed theoretically [14-16]. 

Here we concentrate on the transverse-momentum and impact-parameter dependence 
of hadron/pion ratios at midrapidity and on the nuclear modification factor in dAu at 
^/s = 200 GeV, using leading-order pQCD. We observe that while pQCD is quite successful 
for charged hadron and pion production at large px in proton-proton (pp) collisions, pro- 
ton production in pp is not well understood using the language of pQCD. In fact, pQCD 
underestimates the p/tt"*" ratio by a factor of 3-10 in pp collisions. This can be attributed 
to the limited information about a nonperturbative ingredient, the fragmentation function 
(FF), in the usual factorization-theorem based treatment of particle production in pQCD 
(see D^ZjQ"^) in eq. (1)). Most of the FF-s are extracted from e+e~ data, where the most 
relevant large-z part of the FF-s is not well constrained [17], and the information on proton 
FF-s is very limited, especially for gluon fragmentation functions. We are aware of only 
one set of proton FF-s from global fitting to date [18]. Compared to proton FF-s, FF-s for 
charged hadron and pion production are better studied, supported by more experimental 
information [19,20]. At the same time, the RHIC data for the p/n ratio are only available for 

< 5 GeV, while /i^/tt ratios are available forp^ to about 10 GeV. Therefore, in this paper 
we focus on the h'^/ir^ ratio. Any insight gained is expected to also help in understanding 
the p/tt ratios. 

In an attempt to provide a more satisfactory description of available data and to mimic 
higher-twist contributions (of order QI/pt, where Qy is the appropriate virtuality [21]), 
many pQCD calculations take direct account of the transverse momentum of partons ( "in- 
trinsic fcr")- This can be accomplished via unintegrated parton distribution functions or, 
more phenomenologically, via a product assumption and a Gaussian transverse momentum 
distribution gCkx) (characterized by the width (k^)) [22,23]. We apply the latter procedure 
in the present paper. Then, for pp collisions, the usual convolution of the standard parton 
distribution functions (PDF-s) /„, transverse momentum distributions gi^ir), partonic cross 
sections da/di, and fragmentation functions (FF) D^jc takes the form 
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where the partonic subprocesses a + 6 — > c+d are summed over, Xa, Xb and Zc are momentum 
fractions, and s, t and u denote the parton level Mandelstam variables. In this paper, we 
use CTEQ5L [24] PDF-s and a set of FF-s [19] which were determined relying also on 
semi-inclusive data. 

Figure 1 shows our results for pp collisions. We apply fixed scales, Q = Q = px/^, ior 
both charged hadron and neutral pion production. The top portions of the Figure compare 
the calculated charged hadron and pion spectra to RHIC data [2,25] for pr > 2 GeV, where 
the pQCD calculation can be used. The data are well described with {k^) = 1.8 ±0.3 GeV^. 
The /i^/7r° ratio is displayed in the bottom panel of Fig. 1. The ratio depends very weakly 
on pt, and its value is close to the peripheral AuAu result 1.6) [10]. The dependence 
of the cross sections on the scale is not weak, limiting the predictive power of leading order 
pQCD. Therefore we do not expect to fit the data in detail and concerntrate on the main 
features. One advantage of calculating ratios is that the hadron-to-pion ratio is not sensitive 
to the scale chosen. The width (k'^) has a weak effect on the hadron to pion ratio as long 
as the same (k^) value is used for charged hadrons and for pions. (Later on in this paper 
we will introduce different values of (k^) for different particle species, but such differences 
should be small without medium enhancement, with small effects on the /i^/7r° ratio in pp 
collisions.) 

It was observed in lower-energy proton-nucleus (pA) collisions that the production of 
hard particles {pt > 2 — 3 GeV) is enhanced more strongly than the naively expected scahng 
with A. This so-called Cronin effect [13] has been explained as initial state scattering or kx 
broadening [26]. The observed Cronin effect is not universal for different particle species. 
In pA collisions at ^/s — 38.8 GeV, a stronger nuclear enhancement was seen in proton 
production than for pions [13]. We argue that the non-universality of the Cronin effect may 
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be explained by the non-universality of the kx smearing. The latter can be understood in 
part by the large difference between pion and proton masses. While the incoming parton 
does not "know" whether it will produce a pion or a proton, it requires a larger s to produce 
a proton than a pion at the same transverse momentum. This correlation connects initial- 
state broadening and final-state fragmentation. Furthermore, as mentioned earlier, kr effects 
partly account for higher- twist contributions of order Ql/p^. Assuming larger Ql leads to 
larger power corrections for protons or charged hadrons than for pions. 

We make this observation the centerpiece of our description of the differences between 
the production of different hadrons in the present study. Our strategy is to keep the other 
features of the model as simple as possible. We therefore write 

(4) AS = ikl) +ChLAB{h), 

{klYAB-{kl)+c^LAB{h) , (2) 

where (A;|.) is the width of the transverse-momentum distribution in pp collisions (here 1.8 
GeV^), and we only distinguish between an average value for charged hadron production, 
{k'^)'^, and a value for pions {k^y, corresponding to the data we wish to consider. In both 
cases, we wrote the /cr-broadening as proportional to the effective length of the medium 
along the path of the parton before the hard collisions, L{b), which depends on the impact 
parameter, b. In pA collisions, Lp^ip) is the average of the effective length, 

L{zA,h)= dz pa{z ,h)/po , (3) 

^— oo 

with [za, b) representing the point of the hard collision and po the average density of the 
target, over the Glauber nuclear thickness function of the nucleus. The quantity (3) inte- 
grated over b is proportional to A^^^. In AB collisions, the effective length for a hard parton 
from nucleus A passing through nucleus B, LAB{b), is the average of L{zb, Sb) over the nu- 
clear thickness function for the collision of nuclei A and B, at the given value of b, which 
is proportional to the probability of hard collisions at a certain vector impact parameter b. 
In this paper we use Woods-Saxon nuclear density profiles, with the parameter values for 
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Au taken from Ref . [28] . The effective lengtli can also be written in terms of the number of 
coUisions suffered by the incoming nucleon, ^{b) [23]. For our schematic average purposes 
we find it more appropriate to use the effective length L{b). 

Here we focus on the difference between the coefficients Ch and c^. If the proton contri- 
bution to charged hadrons varied with impact parameter, then the composition-dependent 
Ch should also be expected to depend on b or pt- However, we will show that the data can 
be fitted well by b and pt independent coefficients. This may be understood in terms of 
a pT-integrated proton contribution to the charged hadron yield, which does not strongly 
depend on b or pr- Thus, we neglect any potential b and pr dependence of Ch- We expect 
c/i to be larger than c^. To calculate hadron and pion cross sections separately, one needs 
to account for additional nuclear effects like e.g. shadowing and the suppression mentioned 
in the introduction. Testing various shadowing parameterizations, we found that the shad- 
owing effects are not important to the /i^/tt" ratio. Similarly, if the suppression factors are 
not too different, we expect an approximate cancellation in cross section ratios. We make 
this assumption in the following calculation. 

Figure 2 shows the hadron-to-neutral-pion ratio for different centralities in AuAu colli- 
sions at = 200 GeV with = 0.13 ± 0.04 GeVVfm and cu = 0.45 ± 0.08 GeVVfm. The 
agreement appears to be satisfactory, and is similarly good for other centrality bins (not 
shown). As we expected Ch > c„. In other words, the effective s (the energy involved in 
the partonic cross section) is larger at the same pr for the average hadrons than for pions. 
Larger s leads to more room for a dynamical intrinsic kr, just like in the Drell-Yan case, 
where the larger Ql of the lepton pair leads to a larger [27]. Therefore we find Ch > c,r 
natural, while we do not have a quantitative understanding of the fact that Ch ~ 3.5c^. 
We noticed that the trend of the h^/n ratio at large px for the 60 — 70% bin is somewhat 
different from the trend of other bins. We hope future data with smaller error bars will 
clarify this point. 

As the coUisions become more central, the effective length increases, and the average 
transverse momentum associated with charged hadrons will broaden more than that of 
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pions. This leads to the enhancement of the charged-hadron-to-pion ratio. More study is 
needed to further understand the origin of the difference between the coefficients Ch and c^^. 

Now, let us apply the above parameters to calculate the RdAu nuclear enhancement 
factors for charged hadron and neutral pion production. It is believed that final state effects 
are much less important in dAu than in AuAu, and initial /cr-broadening may play a key 
role. We carry out the calculation with both unmodified PDF-s [24] and nuclear PDF-s 
which incorporate "shadowing". Here we use the EKS parameterization [29]. The results 
are displayed in Fig. 3, together with the recently released PHENIX data [8]. One important 
feature at 2 GeV < < 6 GeV is that the RdAu nuclear modification factors for charge 
hadrons (solid line: without shadowing; dashed: with EKS shadowing) and neutral pions 
(dotted: without shadowing; dot-dashed: with EKS shadowing) are different. The EKS 
shadowing increases RdAu somewhat, since we are in the anti-shadowing region. The nuclear 
modification factors for pions are generally much smaller than the ones for h^, in agreement 
with the data [8]. This is because, as observed e.g. in Ref. [14], the proton- proton (k^) 
determines the position of the Cronin peak, while the coefficient c in eq. (2) regulates the 
height of the peak. The non-universal fc-r-broadening may further contribute to the different 
behavior of the nuclear modification factor for and for pions in AuAu collisions. It will 
however be interesting to carry out a less schematic calculation to understand the different 
shapes of Rab for charged hadrons and neutral pions, respectively. In addition, data for 
the nuclear modification factors are becoming available species by species. A less schematic 
calculation should also address them to distinguish mass effects and effects of particle species 
[30]. 

In summary, we explored the consequences of different Cronin enhancement coefficients 
in the widths of the transverse-momentum distributions of different secondaries in this pa- 
per. We argued for this non-universality on the basis of lower-energy observations and on 
theoretical grounds, related to the higher-twist structure of pQCD calculations. It was found 
that the main features of the transverse-momentum and impact-parameter dependence of 
hadron/pion ratios at ^/s = 200 GeV can be obtained in a simplistic model if the non- 
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universality is incorporated. The different nuclear modification factors recently measured in 
dAu collisions are also displayed by the model. Enhancement of the charged-hadron/pion 
ratio in dAu compared to pp collisions is a direct prediction of this model, which could 
be tested in the near future. The understanding of further intriguing features of nuclear 
modification factors in AuAu collisions is left for future work. 

We are grateful to D. D'Enterria, J. Jia, D. Keane, and J. W. Qiu for stimulating 
discussions. X. F. Zhang acknowledges the hospitality of the Institute for Nuclear Theory, 
where part of this work was carried out. This work was partially supported by the U.S. 
DOE under DE-FG02-86ER-40251. 
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FIGURES 




FIG. 1. Invariant cross section of charged hadron (top) and pion (middle) production in pp 
collisions at ^ = 200 GeV. The data points are from STAR [2] and PHENIX [25]; the solid lines 
represent the leading-order pQCD calculation with {k"^) = 1.8 GeV^. Bottom: hadron/pion ratio 
calculated with (A;^) = 1.8 GeV^ for both, pion and charged hadron production. 
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FIG. 2. Hadron-to-pion ratios at different centralities. From top to bottom: 0-10%; 20-30% ... 
(solid lines); data are from PHENIX(normalization errors are not included) [10]. 
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FIG. 3. The nuclear modification factor RdAu for h (solid line: without shadowing; dashed: 
with EKS shadowing) and vr'^ (dotted: without shadowing; dot-dashed: with EKS shadowing. Data 
points are from PHENIX dAu data [8] . 
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